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Abstract 



We test the convergence of the QCD exponential sum rules by including PT corrections to order a], and the NP contributions up to 

dimension D - % condensates. Then, using the ratio of exponential sum rules where the QCD PT series is more convergent, we 

7— I study the correlation between the gluon condensates (ffjC^) and {g^fabcG^}. From charmonium systems and using the charm quark 

^ mass as input, we deduce: {g^fabcG^) - (8.2 + 1.0) GeV^ x (ffjG^) corresponding to (a^G^) - (7.5 + 2.0) x 10^^ GeV"*. Using these 

results for the bottomium systems, we obtain: mi,(mh) - 4212(32) MeV, which is slightly higher but consistent within the errrors 

with the ones from 2^-inoments and their ratios nihimh) - 4172(12) MeV. We are tempted to consider as a final result from the sum 

O rules approaches, the average mi,(mi,) = 4177(1 1) MeV of the two previous determinations. 
<D \ 
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1. Introduction and motivations 



' _( ' In recent letters llll|2|], we have used finite n and -moments 



- r 

J 4m 



R(t, m]) 

(f + e2)«+i ' 



(1) 



for extracting the values of the gluon condensates and precise 
[ — values of the running c, b quark masses in the MS -scheme. The 
results for the condensates are: 



in 
o 



X 



{a,G^) = (0.07 ± 0.01) GeV* , 
ig'fahcG^) - (8.8±5.5)GeV2x<a,G2), 



(2) 



— I land for the quark masses (average of the two consistent deter- 
^ minations in ^ and iQl): 

mdmc) = 1262(17) MeV , 

mbimh) = 4172(12) MeV. (3) 



Following the same lines, we study in this letter the QCD- 
exponential moments: 



£p (t) = dt t'' e-"R,.At, ml) for p^Q, 

JAml 



and its ratio: 



7 p 

KoiT) = -— Log[i:o] = , 
dr Xo 



(4) 



(5) 



when higher order PT and NP terms are included into their 
QCD expressions; is the ratio of the total cross-section 



of cr(e+e" hadrons) over cr(e+e" — > fj.'^fi^), which is normal- 
ized asymptotically to one; t is the sum rule variable. 
However, working with higher values of p does not help due 
to the relative increase of the higher energy contributions to the 
spectral integral which are not under a good control and which 
decreases the accuracy of the predictions. Then, we shall limit 
ourselves to the case p = 0, 1 where Xo can be derived naturally 
from the moments M„ [O^) by taking the limit: 



oo and ^ oo 



but keeping — - = t finite . (6) 
Q 
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Therefore, by working with Xo and Xi, one can explore a new 
region of energy which is not reached with M„ {Q^^ at finite 

and relatively small n and values like e.g. the M„ {Q^ - O) 
moments. One can also notice that in the case of the exponen- 
tial sum rules, and n are correlated via the sum rule variable 
T, which is not the case of M„ (Q^) where they are completely 
free and uncorrected. In this sense, the two approaches are dif- 
ferent but complementary and one needs to study both methods 
for a complete test of the sum rule approaches. 
Originally named Borel sum rules by SVZ because of the 
appearance of a factorial suppression factor in the non- 
perturbative condensate contributions, it has been shown by 
that the PT radiative corrections satisfy instead the properties of 
an inverse Laplace sum rule. Exponential sum rules have been 
used successfully by SVZ for light quark systems f^. Espe- 
cially, the ratio of moments has been largely used in the litera- 
ture as it appears to be a useful tool for extracting the masses of 
lowest hadrons Jsl-Ql- 

In the case of heavy quarks, Bell-Bertlmann have used the ex- 
ponential sum rules (so-called magic moments) in a series of 
papers for their relativistic and non-relativistic versions l^sl- liol] . 
In particular, they have studied its quantum mechanical inter- 
pretation using the model of harmonic oscillator where they 
found that the sum rule variable t is related to the imaginary 
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time while the optimal value of the results can be extracted at 
the minimum or inflexion point (stability point) of the theoret- 
ical curves. Using this method, it has been emphasized by 
[iQil that the value of the gluon condensate {usG^) obtained by 
SVZ has been underestimated by about a factor 2. A result 
which has been supported by forthcoming papers in difl'erent 



channels Illl4l4l. 



Unfortunately, the existing diff'erent applications of the expo- 
nential Borel/Laplace sum rule for heavy quarkonia systems 
[14I1 and for heavy exotic states (see e.g. ifisll ) suff'er 
from the ambiguous definition of the heavy quark masses used 
in the analysis because, in these applications, the QCD series 
are only known to lowest order or at most to order Us- This 
feature does not permit to extract accurate predictions of the 
quark and hadron masses from exponential sum rules, with the 
exception of working with double ratio of sum rules (DRSR) 
introduced in III6II for extracting mass-splittings and ratio of 
couplings and of form factors despite the poor knowledge of 
the QCD corrections in these channels because the systematics 
of the method and QCD corrections tend to disappear ifisl [PtI - 

H. 

In this letter, we shall improve the previous old extractions 
done 20-30 years ago of the PT (heavy quark masses) and NPT 
(gluon condensates) QCD parameters liMElBby includ- 
ing into the exponential sum rules, pQCD corrections up to or- 
der a] and non perturbative QCD condensates up to dimension 
D - % condensates. Unlike Refs. we shall not also 

use the input value of the {g^fabcG^) deduced from the Dilute 
gas Instanton estimate but leaves it as a free parameter in the 
analysis. In particular, it is important to check these new higher 
order eff'ects on the estimate of the gluon condensate {QsG^) 
and of the quark masses done previously by Bell-Bertlmann Jsl- 
[loll . which we shall compare with some new results obtained in 
|[lll2l using finite and n-moments A1„ (G^)- Indeed, a reli- 
able control of these QCD parameters are mandatory for further 
applications of the SVZ sum rules. 

By parametrizing the high-energy part of the spectral function 
from a threshold tc by a "QCD continuum" which comes from 
the discontinuity of the Feynman diagrams contributions to the 
two-point function, one can transform the previous sum rules in 
Eqs. © and Q into Finite Energy Sum Rules (FESR): 

4 W = r ^'''Re*e-(t, ml) , l^^iT) ^^^Ml (7) 

J Ami. 



where Mr is the lowest resonance mass. 



-'-(1 



2. Expressions of the sum rules 

We parametrize the spectral function using the Jj^ family 
masses and widths for the c-quark channel and the ones of the 
T family in the /7-quark channel. These inputs are given in Ta- 
bles 1 and 3 of We add to these resonance contributions 
the one of the QCD continuum above respectively the threshold 
Vf7 = (4.6 ±0.1) and (11. 098 ± 0.079) GeV, which is described 
by the asymptotic pQCD expression in the massless quark limit 



in the MS' -scheme, which reads for «/ = 4 flavours Q: 

Re*e-\cont ^l+a, + \ .5a] - \2.Qla] , (8) 
where 0: 

as = as(mQ)/7T . (9) 

The perturbative QCD expression is deduced from the well- 
known spectral function to order Os within the on-shell renor- 
malization scheme |25, 2^, which we have transformed to the 
MS -scheme using the relation between the on-shell and MS - 
mass truncated at the same order of PT series ll27ll : 



Mq = niQimQ) ( 1 + -flj 



(10) 



We add to this expression the P T ap proximate spectral function 
to order in the M5 -scheme ll28ll and we add to it the one of 
al term deduced from the expression given in fl^. 
The D = 4 gluon condensate contribution comes from SVZ ^ 
to LO and from ||30|] to NLO. 



The D - 6 and 8 contributions come from Oil 13211 . Their an- 
alytic contributions to the the ratio of exponentiel moments are 
given by lilOtl . We shall use the modified factorization of the 
D = 8 proposed by based on 1 jN expansion with a value 
of about 1/2 of the one from a naive factorization ||32[ 13411 which 
is favoured by the analysis of the Q^-moments in 

3. Tests of the convergence and choice of the sum rules 

In the following, we test the convergence of the charmonium 
sum rules which will permit us to select the quantity which can 
provide reliable results. 

Given the value of mc('We) in Eq. O), one can deduce the nu- 
merical PT series normalized to the LO contribution at the sub- 
traction point v = and taking e.g. t - 1.35 GeV where the 
ratio of sum rules stabilizes in r (see Fig. |2|: 



ZQ{-ml)\pT ^ 1 + 22.6fl, + 85.2fl2 + 152.8a^ 



(11) 



where one can notice that the PT corrections are large. For the 
ratio of moments normalized to the lowest order expression, we 
obtain within the same approximation and for r = 1.35 GeV"^: 



noiJnJpT ^ 1 - 0.12flj - 0.17fl^ + O.Slai 



(12) 



where in this case, the PT corrections are much smaller and the 
series converge quite well. We illustrate the previous discus- 
sions by showing the r-behaviours of the PT moments and of 
their ratio in Fig. [1] We also show the experimental part of the 
moments (continuous-black curve). We have used the central 
values of the input QCD parameters in Eqs. ^ and (|3]l and the 
value of Us deduced from the one from r-decays ll35ll36ll : 



asQuc) = 0.408(14) . 



(13) 



'original papers are in Refs. 317 to 321 of the book (sjl. 
^In the following, we choose the subtraction point v = THq, where all quan- 
tities will be evaluated. 
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Figure 1: a) Behaviour of the moment Xq(w1c) in GeV^ versus r in GeV-2 at 
different orders of perturbation theory: LO: green (thick continuous); NLO: red 
(long dashed); NNLO: blue (short dashed); NNNLO: red-wine (dot-dashed). 
The black (continuous) curve comes from the central values of the data includ- 
ing the QCD continuum, b) The same as in a) but for the ratio of moments 



At this value of r, the non-perturbative contributions to the ratio 
of moments normahzed to the PT LO term read: 



(14) 



which also indicates a quite good convergence of the OPE de- 
spite this relatively large value of t. 

Therefore, in the following, we shall definitely work with the 
ratio of moments for extracting the QCD parameters. 



4. {g fabcG ) versus {asG ) from charmonium 

In this section, our aim is to improve the estimate of the triple 
gluon condensate (g^fahcG^) from the ratio of the exponen- 
tial sum rules using its correlation with the gluon condensate 
{QsG') by giving as input the precise value of the charm quark 
mass rhcinic) in Eq. (O) obtained in IJl, 01 from moment sum 
rules. 

We show in Fig. |2l the r-behaviour of -y/^(m^) by using the 
central values of the QCD parameters in Eq. (I3) and the best 
central value of the instanton radius pc - 1.89 GeV ' obtained 




— — NLO D4 



■ Exp 

■ Expa 
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GeV"-. The two black curves almost degenerate are the data including error 
bars. The green (thick continuous) is the PT contribution including the Z) = 4 
condensate to LO. The long dashed (red) curve is the contribution including the 
Qs correction to the D = 4 contribution. The short dashed (blue) curve is the 
QCD expression including the D = 6 condensate and the dot-dashed (red- wine) 
is the QCD expression including the Z) = 8 contribution. 



in Eq. ( [TtI i by adjusting the minimum of -y/^(m^) to the ex- 
perimental data. One can notice that at the minimum, one has: 



(15) 



indicating that the higher state contributions to the ratio of mo- 
ments are neglible at this value of r. 

We show in Fig. |3]the correlation between pc and {asG^} from 
which we can find a slight stability (inflexion point) for the 
range: 

-2 , 



{asG-} = (7.5 ± 2.0) x 10"^ GeV"* 



(16) 



which is comparable with the one in Eq. ^ obtained from the 
2^-moments though less accurate here. In this region, we can 
deduce in GeV 'Pf 



p, = 1.89(10)c2 (1)„, (5),(1), = 1.89 + 0.11 



(17) 



where the error comes respectively from {a^G^), as, the choice 
of subtraction point v, the localisation of the minimum : t ^ 
(1.25 ~ 1.45) GeV"^. We have estimated the effect of the sub- 
traction point V by varying it from mdmc) to Mj and using the 

Q'i(v) . V 



varying it 
I. nil): 



substitution (see e.g. 

asiific) ai(v)xll-y8i 



log — 

Iflr 



(18) 



where/?! - -(1/2)(11 - 2n//3) for n/-flavours. Theerrordue 
to the a"-"* is obtained by taking the + sign in the coefficient of 
a\ and is negligible. The error due to the data is also negligible. 
The previous value of pc in Eq. (fTTl l leads to: 



ig'fabcG') = (8.2 ± 1.0) GeV^ x {a^G') , 



(19) 



where the relation from the Dilute Gas Instanton approxima- 
tion: 



4 12;r 



(20) 



^If we have used the value of {a^G^) in Eq. (2) from the moments, 
would have obtained a more precise value of p^. 



3 



has been used. We consider this resuh as an improvement of 



the one in Eg. (I2) 
the past HM 



from the 2 -moments 



Notice that in 

1411 . the contribution of {g^fatcG^) has been 
neghgible because of the (a priori) use of the instanton radius 
Pc - 5 GeV ' from the dilute gas approximation fl, which is 
not favoured by the previous result. 
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Figure 3: Variation of pc in GeV ' versus (0,0^) in GeV* from the ratio of 
moments K^ijn^) 



5. jMi(jMi) from bottomium 

Due to the smaller values of asinit,) and of t at which the ratio 
of sum rules stabilizes (see Fig. |4]i, the convergence of the QCD 
PT and NP series is faster here than for the case of charmonium. 
More explicitly, for r = 0.2 GeV"^, the PT expression of the 
ratio of sum rules reads numerically: 

%)(ml)\pT ^ 1 - 0.14flj - 0.35fl? - 0.16al , (21) 

where we have used the QCD inputs in Eqs. (|2]l, Q and (fTTl l. 
We shall use the value of as deduced from the one from r-decay 



Osimt) = 0.219(4). 



(22) 



Using {usG^) from Eq. (|2]i, the improved value of the triple 
gluon condensate in Eq. (fTTT i from exponential charmonium 
sum rule and the modified factorized value of the D = 8 di- 
mension condensates, we reconsider here the determination of 
mhinib) from bottomium by adjusting the minimum (inflexion 

point) of the QCD value of yj^(ml) with the one from the 
data. This is achieved by the value of the running mass in GeV: 

mi(mfc) =4212(18),(0.2)„,(25)v(1.5)g2+c'+c^ (8)co»/, (23) 

corresponding to r = (0.20 ± 0.06) GeV"^ where the inflexion 
point of the theoretical curve meets the experimental one (see 
Fig. H. 

We notice that the sum of the NP condensate contributions is 
about 1 .9 per mil of nih which is quite small. 
The error due to the a"-'* is obtained by taking the + sign in the 
coefficient of al . The one due to v is induced by moving it in the 
range (nii, ±1) GeV. The one due to is the deviation from the 
modified to the usual factorization. The one from Coulombic 
corrections has been estimated like in |2] from the resummed 



expression of the spectral function [i37l] and adding to it the PT 
expression up to order a, computed in OSll which contains the 
familiar (1 - ACfOs) factor due to the quarkonium annihilation 
through a single (transverse) virtual gluon and corrections to 
order v and logv for small v. Then, we deduce the final value: 



mh(mt) ^ 4212(32) MeV . 



(24) 



which is slightly higher and less accurate than the one from the 
2^-moments given in Eq. (OJ ([llSl but the two results are con- 
sistent within the errors. 

By comparing the two results from the 2^-rnoments and ex- 
ponential sum rules, we are tempted to take the average of the 
results in Eqs. ^ and (l24l i as a final result from the sum rules 
approaches: 



mi(mfc) =4177(11) MeV . 



(25) 
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Figure 4: Behaviour of the ratio of moments ^^("Jj) in GeV 

GeV"~ and for rni,(mi,) = 4212 MeV. The same legend as in Fig. [2] for the 
QCD contributions. The black continuous (rep. short dashed) curves are the 
experimental contribution including (resp. without) the QCD continuum. 



6. Conclusions and outlook 

We have reconsidered the QCD-exponential sum rules by in- 
cluding higher order PT corrections to order a] and non pertur- 
bative gluon condensate contributions up to D - 8. We have 
also included the radiative corrections to the D = 4 contribu- 
tion. By studying the convergence of the PT series, we have 
concluded that the ratio of sum rules Ko in Eq. ^ is more con- 
vergent than the moment Xo in Eq. (|4|. 

From charmonium ratio of sum rules ^^(ml), we have deduced 
the ratio between the gluon condensates {(XsG^) and {g^fatcG^} 
with the results in Eqs. (fTTl l and (fT9l l, where we have used 
as inputs the values of mdmc) from Q^-moments in Eq. (|3]l 
Ol and the gluon condensate (aXf') from Q^-moments (HI] and 
from some other sources lll2[|l3ll . We have also used the mod- 
ified factorization for the Z) = 8 condensates proposed by [33i], 
which was favoured in the analysis of Q^-moments The 
present result improves earlier ones obtained in Eq. (13) from 



2 -moments 1 1]. 

Using the previous results to the bottomium systems, we ob- 
tain, from '??Q(m^), the running /j-quark mass in Eq. ( l24l l from 



the ratio of exponential sum rules. The result is slightly higher 
than the one from Q^-moments in Eq. (O and less accurate but 
agrees with it within the errors. We are tempted to consider as 
a final result of Jn[,{mb) from the sum rules approach, the aver- 
age given in Eq. dZSl l from the ratio of Q^-moments and that 
of exponential sum rules. The previous result improves older 
findings for the /7-quark mass from bottomium systems by the 
author ll24ll . by Bell-Bertlmann in a series of papers llsl-liot and 
byH. 
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